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Abstract. We study effects of random fluid motion on a system in a self-organized 
critical state. The latter is described by the continuous stochastic model, proposed 
by Hwa and Kardar [ Phys. Rev. Lett. 62 : 1813 (1989)]. The advecting velocity 
field is Gaussian, not correlated in time, with the pair correlation function of the form 
oc S(t — t')/k d ^ 1+ ^, where k± = |kj_| and ku is the component of the wave vector, 
perpendicular to a certain preferred direction - the d-dimensional generalization of the 
ensemble introduced by Avellaneda and Majda [ Commun. Math. Phys. 131 : 381 
(1990)]. Using the field theoretic renormalization group we show that, depending on 
the relation between the exponent £ and the spatial dimension d , the system reveals 
different types of large-scale, long-time scaling behaviour, associated with the three 
possible fixed points of the renormalization group equations. They correspond to 
ordinary diffusion, to passively advected scalar field (the nonlinearity of the Hwa- 
Kardar model is irrelevant) and to the “pure” Hwa-Kardar model (the advection 
is irrelevant). For the special choice £ = 2(4 — d )/3 both the nonlinearity and the 
advection are important. The corresponding critical exponents are found exactly for 
all these cases. 
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1. Introduction 

Numerous physical systems reveal self-similar (scaling) behaviour over the extended 
ranges of spatial or temporal scales with highly universal exponents. The most 
prominent example is provided by the critical behaviour of various equilibrium systems 
near their second-order phase transition points; see, e.g., CD El and references therein. 
An essentially different example is given by the phenomenon of self-organized criticality 
(SOC), typical of open nonequilibrium systems with dissipative transport; see [3j and 
references therein. In contrast to equilibrium systems, they do not include a tuning 
parameter (like the temperature for the second-order transitions) and evolve to the 
critical state as a result of their intrinsic dynamics. Self-organized critical states are 
believed to be ubiquitous in the Nature and are observed also in biological, ecological 
and social systems [3]. 

Investigation of the effects of various kinds of deterministic or chaotic flows on 
the behaviour of critical systems (like liquid crystals or binary mixtures near their 
consolution points) has shown that the flow can destroy the usual critical behaviour: it 
can change to the mean-field behaviour or, under some conditions, to a more complex 
behaviour described by new non-equilibrium universality classes 0-ra- 

In this paper, we study effects of turbulent advection on a system in a self-organized 
critical state by means of the field theoretic renormalization group (RG). As a rule, SOC 
is studied on the base of discrete models with discrete time steps (like cellular automata). 
We employ the continuous model proposed in nu, described by the stochastic equation 
for a smoothed (coarse-grained) height field h(x) = h(t, x) of a certain profile (sand pile 
or a landscape). Some modifications of the model nn were proposed in [12], [13] and 
the effects of a flow on the erosion of landscapes were discussed in [ISj, but the velocity 
field as an independent dynamical variable was not introduced there. 

In this paper, we explicitly introduce coupling with the velocity field. It is modelled 
by the strongly anisotropic Gaussian ensemble, with vanishing correlation time and 
prescribed power-like pair correlation function, the d-dimensional generalization of 
the ensemble introduced and studied in [13] and, at the same time, the anisotropic 
modification of the popular Kraichnan’s rapid-change model; see [15] for the review and 
the references. It is natural here to consider anisotropic flow, because the original model 
m already involves intrinsic anisotropy, related to the overall tilt of the landscape or 
the sand pile surface. 

2. The model. Field theoretic formulation and renormalization 

Let n be a unit constant vector that determines a certain preferred direction. Then 
any vector can be decomposed into the components perpendicular and parallel to n, for 
example, x = xj_ + nxy with xi • n = 0. Let <9* = d/dxi with i = 1... d be the derivative 
in the full d-dimensional x space, <9j_ = d/dx±i with i — 1... d — 1 is the derivative in 
the subspace orthogonal to n, and <9y = n • d. 
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The stochastic differential equation for the height held h{x ) = h(t, x) is taken in 
the form HU: 

d t h = z/_l 0 d]_h + u\\ 0 dpi - d\\h 2 /2 + /, (2.1) 

where d t = d/dt : u± 0 and up are the bare diffusivity coefficients and f(x ) is a Gaussian 
random noise with zero mean and the pair correlation function 

(/M/M) = 2 D 0 S(t - f) «W(x - x'). (2,2) 

Coupling with the velocity held vf x) is introduced by the replacement d t —>• V t = 
dt + lydi, where V* is the Lagrangean (Galilean covariant) derivative. 

The velocity held will be taken in the form v = rry(t,xj_), where r>(t,xj_) is a 
scalar function independent of x\\. Then the incompressibility condition is automatically 
satished: 

diVi = d\\v(t,x ± ) = 0. (2.3) 

For u(t,xj_) we assume a Gaussian distribution with zero mean and the pair correlation 
function 

[' rllr 

(v(t,x ± )v(t',x.' ± )) = 5(t-t') J expjik- (x-x)}A,(fc) = 

= 5(t-t') J ^a-i exp H k -L ' (x± - x'_l)} D v (k ± ), (2.4) 

with k± = |kj_| and the scalar coefficient functions of the form 

D v {k) = 2n6{k\\) D v (k ± ), D v (k±) = B 0 kf d+1 ^. (2.5) 

Here Bq > 0 is a constant amplitude factor and f is an arbitrary exponent, which 
(along with the conventional e = 4 — d) will play the role of a formal RG expansion 
parameter. The infrared (IR) regularization in (12.4p is provided by the cutoff k± > m. 
[Its precise form is unimportant; sharp cutoff is merely the most convenient choice from 
the calculational viewpoints.] 

According to the general statement (see, e.g, the monographs m ei and the 
references therein), our stochastic problem is equivalent to the field theoretic model 
of the extended set of fields 4> = {h', h, v} with action functional 

5(4*) = h'D 0 h' + h!{—d t h — vd\\h + u± 0 d 2 ± h + upd 2 \h — d\\h 2 /2} + . (2.6) 

All the required integrations over x = { t , x} and summations over the vector indices 
are implied. The last term in (I2.6[) corresponds to the Gaussian averaging over v with 
correlator (j2.4j) and has the form 

^ J dt J dx.±dx' ± v(t, x_ l )5“ 1 (xj_ — x' ± )v(t, x'j_), (2.7) 

where 

5„- 1 (r ± )<xD 0 -‘rf-(2.8) 


is the kernel of the inverse linear operation D v 1 for the correlation function D v in (12. 5 p . 








Effects of random environment on a self-organized critical system 


4 


This allows one to apply the field theoretic renormalization theory and 
renormalization group to our stochastic problem. The model (12.6[) corresponds to a 
standard Feynman diagrammatic technique with three bare propagators: (vv)o, given 
by (12.4ft . 112.5]) . and the propagators of the scalar fields (in the frequency-momentum 
representation): 

{hh') 0 = {h'h}* 0 = {-ioj + e(k)}- 1 , (hh) o = 2D 0 {uj 2 + e 2 {k )}~ 1 , (2.9) 

where e{k) = v^kf + "lo k l- The propagator ( h'h') 0 vanishes identically for any field 
theory of the type (12.61) . The model also involves two types of vertices corresponding to 
the interaction terms — h'd\\h 2 /2 and — h'(vdn)h. The corresponding coupling constants 
go and wq are defined by the relations 

D 0 = <4%, B 0 = w 0 v\\o, ( 2 . 10 ) 

so that by dimension g 0 ~ £~ e and Wo ~ £~^, where is £ has the order of the smallest 
length scale in our problem. 

The analysis of canonical dimensions (see, e.g., m Ei) along with symmetry 
considerations shows that our model is multiplicatively renormalizable with the only 
independent renormalization constant Z v ^\ all the renormalization constants are equal 
to 1. In the leading (one-loop) approximation Z u has the form =1 — ag/e — bw/f 
(the use of the minimal subtraction scheme is implied). Here g and w are renormalized 
analogues of the bare parameters in ( 12 . 101 ) . which play the role of the coupling constants 
(dimensionless expansion parameters) in renormalized perturbation theory, and a, b 
are numerical coefficients. Their precise values are unimportant (they can simply be 
absorbed by the redefinition of g and w), but it is important for the following that they 
are positive: a, b > 0 . 

To proceed, we have to find the RG functions: the anomalous dimensions 
and the (3 functions, which appear as coefficients in the differential RG equations; 
detailed discussion can be found in mq. The anomalous dimension 7 p for a given 
renormalization constant Zp is defined as 

7 f = V M In Zp for any quantity F, (2.11) 

where V^ L is the differential operation /ad^ for fixed eo, the latter being the full set of bare 
parameters eo = {^,^± 0 ,^ 0 ,9o}- The f3 functions for the two dimensionless couplings 
g and w are 

P g = D^g = g [-e - 7 5 ], = V^w = w [-£ - y w ]. (2.12) 

Taking into account the relations between the bare parameters (12.10]) . absence of 
renormalization of D 0l B 0 and the explicit form of the constant Z v , we obtain: 

A? = 9[-£ + 37 „n/ 2 ], /3 w = w[-f + 7 „ n ], 7q, = ag + bw + ... (2.13) 

where the cllypsis stands for the higher-order corrections in g and w. 
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3. Fixed points and scaling regimes 

It is well known that possible scaling regimes of a renormalizable model are associated 
with IR attractive fixed points of the corresponding RG equations. In the present case, 
the coordinates g*, tc* of the fixed points are found from the equations 

/3 g (g*,w*) = 0, /3 w (g*, wf) = 0, (3.1) 

with the fd functions given in (12.131) . The type of a fixed point is determined by the 
matrix = {f= dfdi/dgj}\ g ,, w ,, where fdi is the full set of the fd functions and 
g 3 = {g, w} is the full set of coupling constants. For IR attractive points the matrix hi 
is positive, that is, the real parts of all its eigenvalues hi, are positive. 

The f functions (12.131) involve the same anomalous dimension 7 ^ and thus satisfy 
the exact relation wfd g — 3g[d w /2 = gw(—e + 3£/2). I 11 turn, this means that for general 
e and £ the equations (13. lj) can be satisfied only if at least one of the coordinates g*, u>* 
vanishes identically. There are three such fixed points: 

(1) Gaussian (free) fixed point with g* = w* = 0 and 7 * = y^g* = w* = 0) = 0. 
The eigenvalues of the matrix hi are equal to its diagonal elements: f 1 g = d g fd g = — e, 
£l w = d w fd w = — £. Thus, the point is attractive for e, £ < 0. 

(2) The point with g* = 0, w* = £/6 and 7 * = £. The latter expression is exact 
(no corrections of order £ 2 and higher) owing to the exact relation between /3 W and 
7^11 in (12.131) . Furthermore, the derivative d w j3 g = 3gd w r ) v j2 vanishes for g* = 0, the 
matrix fl is block triangular and its eigenvalues coincide with the diagonal elements: 
Q g = — £ + 3£/2 (exact result) and tt w = £ (with possible corrections of order £ 2 and 
higher). We conclude that this point is attractive for £ > 2e/3, £ > 0. 

From the physics viewpoints, in this critical regime the nonlinear term in the 
equation (12. lj) does not affect the leading terms of the IR asymptotic behaviour (it 
is IR irrelevant in the sense of Wilson). 

(3) The point with w* = 0, g* = e/a and 7 * = 2e/3. The latter expression is 
exact (no corrections of order e 2 and higher) owing to the exact relation between fd g 
and 7 ^n hr (j2.13j) . Furthermore, the derivative d g /3 w = wd g -) v ^ vanishes for w* = 0, 
the matrix hi is block triangular and its eigenvalues again coincide with the diagonal 
elements: £l w = —£ + 2e/3 (exact result) and Q w = e (with corrections of order e 2 and 
higher). We conclude that this point is attractive for £ < 2e/3, e > 0. In this critical 
regime, the advection is irrelevant and the IR behaviour of the model coincides with 
that of the original Hwa-Kardar model. 

It is also worth noting, that for all fixed points (l)-(3) in the regions of their IR 
attraction, their coordinates lie in the physical region g* > 0 , w* > 0 . 

For the special choice £ = 2e/3 the functions (I2.13j) become proportional, and 
equations (13. ip only impose the restriction ag* + = £; that is, we have a line of fixed 

points in the g-w plane. However, the exact value 7 * = £ = 2e/3 is the same for all 
points on that line. The direct calculation of the elements of the matrix shows that 
one of the eigenvalues vanishes (this exact result is a consequence of the degeneracy of 
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Figure 1. Regions of attraction of the fixed points in the model (12.61) . It should be 
noted that all the boundaries of those regions are straight rays; there are neither gaps 
nor overlaps between the different regions. These results are exact (not affected by the 
higher orders £ 2 , e 2 and so on.) 


the fixed point) while the other, equal to 3ag*/2 + bw*, is positive in the physical region 
> 0, as well as the exponent £. 

In figure [1] we show the regions of IR attraction for all the fixed points in the 
plane, that is, the regions where the eigenvalues of the matrix are both positive. 

4. Critical scaling and critical dimensions 

Existence of IR attractive fixed points of the RG equations leads to appearance of the 
scaling behaviour of the Green functions in the IR range. In dynamical models with two 
momentum scales, like our present model, the critical dimension Ap of an IR relevant 
quantity F (a field or a parameter) is given by the relation 

A p — dp + d|'A|| + d°f A^ + 7 p. (4-1) 

Here dp are the canonical dimensions (see, e.g., m), 'yp is the fixed-point value of the 
corresponding anomalous dimension, 

Ay = l + 7 *,/2, A £J = 2- 7 : x (4.2) 

are the critical dimensions of the momentum in the direction parallel to n and of 
the frequency, respectively, with the normalization Ay = 1 for the momentum in 
perpendicular subspace. Detailed derivation and discussion of the relations of the type 
(14. ip . (I4.2[) can be found, e.g., in ref. |9]. 

In our model these relations simplify owing to the absence of renormalization of all 
the fields (so that 7 ^ = 0) and of the parameter up 0 (so that 7 * x = 0); see sec. [2] What 
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is more, 7 * is known exactly for all the fixed points. This gives the following exact 
expressions: 

A h = 1 , A h , = d = 4 - e, A w = 2 , A„ = 1 (4.3) 

for the Gaussian fixed point (1), 

A h = l-f/2, A h , = d + f = A-e + f 1 A u = 2 , A,, = 1 + £/2 (4.4) 
for the hxed point ( 2 ) and the line of the fixed points, and 

A ft = 1 -e/3, Aft/ = d + 2e/3 = 4 — e/3, A w = 2, Ay = 1 + e/3(4.5) 
for the point (3). 

In particular, for the pair correlation function of the main held this gives: 

(h(t, x) h(0, 0 )) ~r" 2Ah E (t/r^fr\\/r^ , (4.6) 

where r± = |xj_| , r\\ = x\\ and T is a certain scaling function of critically dimensionless 
arguments. 

It remains to note that for the hxed point (3) our results agree, up to the notation, 
with those obtained in mi (one has to identify z = A^/Ay, £ = l/Ay, x — — Aft/Ay). 

5. Conclusion 

We studied ehects of fluid motion, including turbulent mixing, on a system in the state 
of self-organized criticality. The latter was modelled by a stochastic differential equation 
for a height held ( 12 . 11 ) . ( 12 . 2 ft . while the velocity of the fluid was described by a random 
Gaussian ensemble (12.41) . Both the equation and the velocity ensemble are strongly 
anisotropic. 

The full model can be reformulated as a renormalizable held theoretic model with a 
single independent renormalization constant. The corresponding RG equations possess 
several IR attractive hxed points, corresponding to different possible types of critical 
behaviour. They correspond to ordinary diffusion, passively advected scalar held and 
to the original SOC model without the mixing. Their regions of IR attraction in the 
plane of the model parameters d and £ and the critical dimensions of the basic helds 
and parameters are found exactly. For the special choice £ = 2(4 — d)/3 an intermediate 
regime arises, where both the nonlinearity and the advection are important. 

In the following, it would be interesting to study other models of SOC and to 
employ more realistic models for the turbulent velocity held, e.g., models with finite 
correlation time and non-Gaussianity. This work is in progress. 
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